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We show how population inversion (PI) occurs in a two-level system (TLS) while measuring its 
charge using a capacitively coupled superconducting single-electron transistor (SSET), biased in the 
subgap regime, where the current through the SSET is carried by different cycles involving tunneling 
of both Cooper pairs and quasiparticles. The PI is directly associated with the resonant nature of 
the Cooper-pair tunneling. We also show how the SSET may strongly relax the TLS, although there 
is negligible current flowing through the SSET, i.e. it is turned off. The calculation of the quantum 
back-action noise is based on a real-time Keldysh approach. We specifically discuss the case of a 
Cooper Pair Box qubit with the SSET capacitively coupled as read-out device. 

PACS numbers: 03.67.Lx 42.50.Lc 73.23.Hk 85.25.Na 



The single-electron transistor (SET) is known as an 
extremely sensitive electrometer jl], |) based on the 
Coulomb blockade ]|, Q. The SET consists of a small 
capacitance island, isolated from the leads by tunnel 
junctions. The measurement is performed by measuring 
the current flowing through the island, which depends 
strongly on the charge induced through the capacitive 
coupling to the charge to be measured. 

With the invention of the RF-SET||, the SET was 
also made fast, and operating frequencies above 100 MHz 
could be reached. Therefore the SET can be used in 
applications ranging from very sensitive charge meters 
and current standards j| in which electrons are counted 
or pumped one by one, to read-out of quantum bitsj^, |7], 
H |9| or to work as photon detectors 

For these sensitive applications the back-action of the 
SET during read-out, arising from the voltage fluctua- 
tions on the SET island, is important since it may corrupt 
the measurement jO, [l^] . 

The SET may also be used in the superconducting 
state (SSET) and indeed such a SSET shows the state-of- 
the-art sensitivity of 3.2- 10" 6 c/Vlfe |J|. This SSET was 
biased around the gap edge, i.e. the threshold voltage for 
sequential quasiparticlc tunneling |fl4f , where the back- 
action may be analyzed in the same way as for sequential 
electron tunneling in the normal state jy], O, . 

There have also been suggestions for using the SSET 
biased in the subgap regime (jl6| |l7]], where the current 
is carried by different sequences of resonant Cooper- 
pair tunneling and quasiparticle tunneling, so called 
Josephson-Quasiparticle (JQP) cyclesjl], [l^] . 

The current noise, determining the shot-noise limited 
sensitivity was recently analyzed for the simplest JQP 
cycle This process consists of a Cooper-pair tun- 
neling across one junction followed by two quasiparticles 
tunneling across the opposite junction. Clerk et al|20| 
analyzed a more complicated process with two Cooper 
pairs involved, called the Double JQP process (DJQP), 
which corresponds to the so called eye-feature in the 



SSET current-voltage characteristics, see Fig. ||b. The 
resonant Cooper-pair tunneling was shown to improve 
sensitivity. 

For the DJQP it was shown that also the back-action 
changes qualitatively and may even induce population 
inversion (PI) in the measured system[p0|. This shows 
that the back-action of the SSET indeed is very non- 
trivial, e.g. simply assigning a noise-temperature to this 
device would demand negative temperatures. 

In this paper we describe the back-action in the whole 
subgap regime. We specifically map out where PI is pos- 
sible, which also includes the simplest JQP process. Fur- 
thermore, we show where the relaxation of the measured 
system may be strong, even though there is negligible 
current flowing through the SSET, i.e. it is turned off. 
We choose the Single Cooper Pair box (SCB)|plfl as our 
specific two-level system (TLS). 

The Model Consider a small superconducting SET 
island coupled via low transparency tunnel barriers to 
two external superconducting leads, and coupled capaci- 
tively to a control gate and to the SCB, see Fig.|[ The 
voltage noise on the SET island is calculated neglect- 
ing the coupling to the SCB. This approach is appro- 
priate in the considered limit of weak SET-box coupling 
(C c « C L ~ C R ~ C b ). 

We follow the outline of Ref. g| and model the SSET 
by the Hamiltonian 

H = H l +H r +Hi+Hc+K+Hj+H t = H +H T , (1) 

where H r = £ fcrf e r kn a\ rj} a krn and Hi = £ in ti n c\ n c ln 
describe the noninteracting quasiparticles in the r € 
{L, R} leads and on the island. Introducing the oper- 
ators fiL/R for the number of charges passed from left 
to right in the L/R junction we write the Coulomb term 
H c = E C {N — n x ) 2 where E c — e 2 /2Cs is the charging 
energy and N = ht — n R is the operator for the excess 
number of charges on the island. Cs ~ Cl + Cr is the 
island capacitance neglecting the small gate and coupling 
capacitances, and n x is the fractional number of electrons 
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FIG. 1: Schematic figure of the superconducting single- 
electron transistor capacitively coupled to a two-level system, 
in this case a Single Cooper Pair box (SCB). 



induced by the gate voltage C g V g /e. The work done by 
the biasing circuit is K — eVd s (CnhL + ClTIr)/^ and 
the Josephson coupling is Hj — — ^ r cos 2<J> r , where 
e r increases n r by one. Then we work in the eigenbasis 
of -f/oj and treat the quasiparticle tunneling Ht as the 
perturbation. We use a real-time diagrammatic Keldysh 
techniquejl^, |2^| to calculate the DC current and the 
spectral density of the SET island voltage fluctuations, 

e 2 f°° - 
Sv{u) = / dre+"- T Tr{p st (t o )N(r)N(0)}, (2) 

which determine the back-action. We consider lowest 
non- vanishing order in Ht, i.e. we neglect co-tunneling 
and higher order effects. Due to the low junction trans- 
parency we also neglect logarithmic renormalization ef- 
fects plf . In this approximation the rates for the 
quasiparticle processes involved are given by the usual 
current- voltage characteristics for tunneling between two 
superconductors p3|, where the effective voltage is given 
by the energy gain in the process. 

Assuming that the SSET is the strongest noise source 
coupled to the two-level system (TLS) the rate for relax- 
ing/exciting the TLS is given byi, M, |25): 



Vt = "72 S v (±AE/h) 



(3) 



where a is a dimensionless coupling constant and AE 
is the energy splitting between the two states in the 



TLS. For the SCB a = %t^±i- To describe pop- 
ulation inversion we calculate the polarization of the 
TLS in the steady state, i.e. AP = P[ LS - P TiS = 
(r-f — r^)/(r-f + r^). The total back-action induced mix- 
ing rate Tf +T| is proportional to the symmetrized noise 
Sy m (AE/h) = Sv(AE/h) + Sv{-AE/Ti). If this rate 
is small, then other noise sources may determine the 
steady-state of the two-level system. To see where the 
back-action is strong we also study Sy m (AE/h). 

Results To be specific we now consider a symmetric 
SSET with the same parameters as in reference |20|, i.e. 
with charging energy equal to half of the superconducting 



energy gap, E c = A s « 25meV, C L = C R , and R T 
1 1 1 « 50kfJ, using the Ambegaokar-Baratoff value for the 
Josephson energy, Ej = Ef ~ 16/xeV. We choose a level 
splitting of the TLS of AE = 0.1A 5 . To describe the 
main features in the current through the SSET and the 
noise in this regime it is enough to take the five lowest- 
energy charge states into account. Numerically we check 
that including 7 states does not change the result. In 
Fig. we show the numerical results for the DC current 
through the SSET, the symmetrized noise S v vm {AE/Ti) 
and the polarization of the TLS AP, as a function of 
SSET gate voltage and driving bias. 

Analytical Results At specific bias points in the sub- 
gap regime, where specific JQP processes are known to 
dominate, one may calculate Sy{uj) analytically. The 
main approximation is to consider all involved quasipar- 
ticle transition rates equal and frequency independent. 
For frequencies not too high compared to the supercon- 
ducting energy gap (fiu) < As) this is a reasonable ap- 
proximation. Keeping track of the density matrix of up to 
5 charge states, and the coherence in two resonant pairs, 
leads to inverting matrices of dimension up to 9 x 9. Us- 
ing Mathematica we obtain analytical results, which are 
presented below at appropriate places. For the DJQP 
cycle, in the proper limit, this procedure reproduces the 
formulas for <SV(w) derived by Clerk et alp(j|]. 

Simple JQP resonances The Coulomb energy of the 
charge state with N extra electrons on the island is 
E N = E C (N — n x ) 2 . Along the lines eVds — E± 2 — E = 
4(1 =F n x )Ec in Fig. || the charge states N = and 
N = ±2 have equal energy, including the work 2(eVd s /2) 
done by the bias voltage across the L/R junction. (Since 
Cl = Cr the voltage divides equally over the junc- 
tions.) Due to the Josephson coupling Cooper pairs tun- 
nel resonantly E3| back and forth across the L/R junc- 
tion. For eVds > 3Ec the resonantly tunneled Cooper 
pair may decay into two quasiparticles tunneling across 
the opposite (R/L) junction, thus completing the sim- 
plest JQP cycle [[jj, which effectively transports two elec- 
trons through the SSET, see Fig. [|a. Approximating the 
two involved quasiparticle rates with a single frequency- 
independent rate T/h and defining the energy gain for 
the Cooper-pair tunneling as S± — eVds — 4(1 =F n x )Ec, 
the population polarization close to a single resonance, 
but not close to both is 
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with maximum max{AP J Q p } = 2AE/\ / T 2 + 4AP 2 lo- 
cated at 5± = ^/(r/2) 2 + AP 2 . We see that for positive 
8± the population is inverted, and that the maximum 
inversion grows with increased TLS level splitting AP. 

Qualitative explanation of the PI Without interaction 
with the TLS the effective Cooper-pair tunneling rate has 
a maximum at zero energy gain (8± = 0)p0|, p6| . In a 
JQP cycle which excites the TLS the effective energy gain 
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FIG. 2: Numerical results for a symmetric SSET with Ec — 
As » 25meV, Rt = Rl « 50kft and Pj = Pj « 16/xeV. 
The level splitting in the two-level system is AP = 0.1 As. 
Top panel - DC current through the SSET. Middle panel - 
Total noise S^ m (AE/h). Lower panel - TLS population po- 
larization AP = Pi — Pq. 
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FIG. 3: a) Two simple JQP processes in parallel. The 
Cooper-pair tunneling in the left/right cycle is resonant along 
the line eVds = 4(1 ± n x )Ec- b) The Double JQP cycle, res- 
onant around eVds = 2Ec and n x — 0.5. 




FIG. 4: Schematic figure of effective tunneling rates for 
Cooper pairs (a-c) and quasiparticles (d) in the SSET, as 
function of the energy gain SE in the process. The vertical 
dashed lines denote the energy gain without transitions in the 
TLS, which in (a-c) is S±. Ff /j, indicate tunneling rates with 
simultaneous excitation/relaxation of the TLS, where the en- 
ergy gain is shifted by AP. 



is lowered to S± — AE, see Fig. ||. For positive S± this 
is closer to resonance (zero) than S± + AE, which is the 
gain for the relaxing cycle. Thus the exciting cycle runs 
faster than the relaxing one, causing a PI of the TLS. As 
for the quasiparticle tunneling rates they monotonously 
increase with increasing energy gain (see Fig. |^d). Thus 
PI will never occur above the threshold voltage, where 
the current is carried only by quasiparticles. There the 
relaxing (pure quasiparticle) cycle will always run faster 
than the exciting one. 

Two JQP cycles in parallel Close to eVds = ^Ec 
and n x = 0, Cooper-pair tunneling across both the left 
and right junctions are resonant simultaneously. Here 
the two JQP cycles run in parallel, see Fig. ||a, shar- 
ing the charge state N = 0. From the numerical re- 
sults we deduce that the maximum population inversion 
is not larger at this crossing compared to the single res- 
onances. Again approximating all involved quasiparticle 
rates with a single T, and now also taking the relevant 
limit Ej <C r, we find along the vertical symmetry line 
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eVds — 4:Eq and n x = 0, that the popula- 



tion is indeed given by Eq. (Eh . Along the horizontal line 
eVds = 4Ec we have <5+ = — <5_ giving AP = 0. 

The Double JQP cycle Around eVds — 2Ec and 
n x = 0.5, Cooper-pair tunneling between the charge 
states N = — 1 and N = 1 across the right junction 
and between N = and N — 2 across the left junction 
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are both resonant processes. This is where the Double 
JQP process |20] occurs, transporting 3e in each cycle, see 
Fig |3|b. The population polarization here is remarkably 
similar to the one around the simple JQP crossing, except 
that here it has a real maximum along the symmetry line 
n x = 0.5. For Ej <C T the maximum is again located at 
6 = v /(r/2) 2 + AE 2 , where now 5 = eV - 2E C , and the 
value is 

max p2 _i_ 8Ai? 2 ' 

Strong Relaxation - No Current For the SSET back- 
action to be important compared to other noise sources 
the total noise Sy m (AE/h) should be large. In the mid- 
dle panel of Fig. || we see that the noise is large where 
there is substantial current, with one exception: The 
noise is also large below the lines eVd s /2 = E±\ — Eq — 
(1 ± 2n x )Ec, where the width of the region is set by the 
TLS level splitting AE. Here a strong relaxation of the 
TLS occurs, although there is negligible current. In this 
bias regime the Cooper-pair tunneling is off resonance, 
but it still slowly takes the SSET from the lowest energy 
charge state N = to higher energy states. Above these 
lines the SSET relaxes through quasiparticle tunneling. 
Thus the SSET mainly occupies its lowest energy charge 
state N = 0. Below these lines this relaxation of the 
SSET is not energetically allowed, and the SSET ends up 
with a population of the excited states N = ±1 of order 
unity. In the region with strong noise the quasiparticle 
transition iV = ±ltoiV = 0is possible through absorp- 
tion of the TLS energy, thus causing strong TLS relax- 
ation although there is negligible current flowing through 
the SSET. In short, the slow off-resonance Cooper-pair 
tunneling pumps the SSET to an excited state which may 
only relax through a simultaneous relaxation of the TLS. 

Conclusion Using a real-time diagrammatic Keldysh 
approach we have calculated the back-action noise of a 
superconducting single-electron transistor, biased in the 
subgap regime. We find that population inversion of the 
measured two-level system is possible when the JQP cycle 
which gives energy to the TLS is closer to resonance than 
the JQP cycle which takes energy from the TLS. We also 
find regions where the relaxation of the TLS is strong, 
although there is negligible current flowing through the 
SSET, see Fig. § 
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FIG. 5: The dark grey indicates for which SSET bias PI of 
the TLS occurs. Thick black lines indicate substantial sub- 
gap current through the SSET. The light shaded areas indi- 
cate where strong TLS relaxation occurs, although there is 
negligible current running through the SSET. 
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